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Abstract 


We explore how nuclear modifications to the nucleon parton distributions 
affect production of high-pj’ hadrons in heavy ion collisions. We perform a 
leading order calculation of the high-py charged pion, kaon and proton spec¬ 
tra using standard fragmentation functions and shadowing parameterizations. 
We also consider alternate models of shadowing. Near midrapidity, shadowing 
is a small effect and cannot explain the large observed suppression of high-p^ 
hadrons. We also consider the isospin difference between protons and nuclei 
and find that it is also a small effect on high-p^ hadron production. 
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Some of the most interesting experimental results to come out of the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Laboratory have involved the observed sup¬ 
pression of high-pr hadron production in collisions of gold nuclei at a center of mass energy, 
VSnn, of 130 and 200 GeV per nucleon. Both the PHENIX [1] and STAR [2] collabora¬ 
tions have published results showing that production of hadrons with transverse momentum 
Pt > 2 GeV/c are suppressed compared to the pp reference spectrum convoluted with the 
number of binary collisions. The suppression is given in terms of the ratio 

daAA/dpT 

TAAdUpp/dpT 

where daAA/dpr and dapp/dpr are the hadron pt distributions in AA and pp collisions 
respectively and Taa is the nuclear overlap function. The product of the inelastic nucleon- 
nucleon cross section with the nuclear overlap function is the number of binary nucleon- 
nucleon collisions in a given impact parameter range. 

The data show that, for p^ > 2 GeV/c, Raa{pt) is much less than 1. For charged 
hadrons (pions, protons and kaons), Raa{pt) ~ 0.4 [1,2], while for 7r°, Raa{pt) ~ 0.3 [1]. It 
appears that, for protons, there is very little suppression, with Raa{pt) ~ 1-0 [3]. A related 
ratio has been formed for central to peripheral collisions which produce low multiplicities, 
such as those in pp interactions. Similar values were found [1,2]. The results at y/S^iN = 130 
GeV and 200 GeV are in agreement with each other. 

One possible explanation for the strong suppression is parton energy loss in the medium 
produced in the collision [4,5]. However, other, more conventional nuclear effects must also 
be considered. Gold nuclei have a different isospin from the proton reference. Isospin is much 
less relevant for comparison of central and peripheral ion collisions. More importantly, the 
parton distributions in nuclei are known to be different from those in bare nucleons [6]. In 
their analysis, the experimenters noted that this difference, referred to as nuclear shadowing, 
might affect Raa but estimated that the change would be small. Here, we give quantitative 
estimates of the effects of nuclear shadowing and isospin on Raa for charged pions, kaons and 
protons separately. The Gronin effect, which broadens the pt distributions in pA relative 
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to pp collisions, tends to increase Raa [4]. To better illustrate the effects of shadowing and 
isospin alone, we do not include the Cronin effect in our calculations. 

We make a leading order (LO) calculation of minijet production to calculate the yield of 
high-pr partons [7]. The pr distribution is [8] 


da^B 


= 2pT 


dxjfi/A{xi, Q‘^)fj/B{xj, Q 


^ Dh/k{zc) da 


where x is the fraction of the hadron momentum carried by the interacting parton and Q 
is the momentum scale of the interaction. The integral over center-of-mass scattering angle 
0 < 9cm < TJ" corresponds to an integral over all rapidities. The 2 —2 minijet cross sections 
are given by da/di. The fragmentation functions, Dh/k{zc), are the probability for the 
production of hadron h from parton k [9]. The parton densities, fi/A{xi,Q‘^), in a nucleus 
are fi/A{xi) = R\{xi){ZAfi/p{xi) + NAfi/n{xi))/A where Za and Na are the proton and 
neutron numbers in nucleus A. We use the MRST LO proton parton distribution functions 
for isolated nucleons [10] and the EKS98 parameterization of the shadowing function R\ 
[11]. We take = Pt- 

In our calculations, we neglect higher-order corrections, using LO parton densities, frag¬ 
mentation functions, and shadowing parameterizations. Since Raa is a ratio, the higher- 
order corrections should largely cancel out. The produced quarks, antiquarks and gluons 
are fragmented into charged pions, kaons and protons using LO fragmentation functions 
£t to e+e“ data [9]. The hnal-state hadrons are assumed to be produced pairwise so that 
TT = (7r+ -|- 7r“)/2, K = (77+ -|- K~)/2, and p = {p + p)/2. Any baryon asymmetry in the 
initial state then has no effect on the baryon composition of the hnal state. The produced 
hadrons follow the parent parton directions. The evolution is modeled using data 
at several different energies. These fragmentation functions were also compared to pp, yp 
and 77 data. After some slight scale modifications [12] they were able to fit all the h~ data. 
However, there are significant uncertainties in fragmentation when the leading hadron takes 
most of the parton momentum [13]. 

Shadowing is assumed to depend only on the parton momentum fraction x, the momen- 
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turn scale Q^, the parton flavor and the nuclear mass number A. The EKS98 parameteri¬ 
zation evolves each parton type separately for 2.25 < < 10^ GeV^. In the momentum 

range relevant to RHIC, R\ for quarks and antiquarks are based on nuclear deep-inelastic 
scattering data. There is very little direct data on the nuclear gluon distribution so that 
gluon shadowing is primarily based on the evolution of the nuclear structure functions. 
The gluon density shows significant antishadowing for 0.1 < x < 0.3 while the antiquark 
densities are shadowed in this region. For 0.3 < x < 0.7, there is significant suppression for 
all partons, the EMC effect, while for x < 0.07, there is also significant suppression. 

We do not include inhomogeneous (spatially varying) modifications [14]. Shadowing 
should be largest near the core of the nucleus and reduced near the surface. In very peripheral 
collisions, shadowing effects should be sufficiently reduced for pp collisions to be a reasonable 
model of the most peripheral events. For the most central collisions, the deviation from 
Raa = 1 differs by less than 1% from than that calculated with homogeneous shadowing. 

Minijets come from quarks, antiquarks and gluons produced in the hard scattering. The 
fragmentation function determines how these partons become final-state hadrons. Figure 1 
shows the percentage of charged pions, kaons and protons produced by quarks, antiquarks 
and gluons as a function of pt at central rapidities, \y\ < 1, for pp interactions at ^/S = 200 
GeV. Pion production is dominated by gluons up to pr ~ 10 GeV/c. At higher px, most of 
the pions come from quarks. For kaons and protons, the crossover between gluon and quark 
dominance occurs at lower px, px ~ 3.5 GeV/c for kaons and ~ 5 GeV/c for protons. At 
large px, ~ 75% of kaons and protons are produced by quarks. At larger rapidities and low 
Px, the relative quark and antiquark contributions increase for kaons and protons. In fact, 
the gluon contribution to kaon production becomes smaller than the quark and antiquark 
contributions. 

Figures 2 and 3 compare the calculated Raa integrated over all rapidities to that re¬ 
stricted to Ip I < 1. The region |p| < 1 roughly matches the experimental acceptances. The 
results shown in Figs. 2 and 3 probe different x regions, as described below. Without frag¬ 
mentation, Xij = {px /where pi and p 2 are the parton rapidities. When 
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yi = y 2 = 0, Xi = Xj = 2pt/^/Snn- For illustration, we assume that the leading hadron of 
the minijet carries half the parent parton momentum, increasing Xi by a factor of two. Then 
the region 2 < px < 10 GeV corresponds to 0.04 < Xi < 0.2. In this x range, the EKS98 
parameterization has valence quark and gluon antishadowing and sea quark shadowing. The 
shadowing modihcations decrease with Thus the strongest effects are at low pt but the 
overall effect is small. Away from y = 0, Xi and Xj are different. One x will decrease into 
the low X shadowing region while the other will increase into the EMC region. There are 
stronger modihcations in both these regions so that shadowing has a bigger effect in the 
broader rapidity region. This is most clearly seen for low pt in Fig. 2 where shadowing 
reduces Raa ~ 30% for all species. At higher pt the effects are smaller and only the proton 
curve comes close to the data. The large suppression observed at RHIC for mesons is not 
seen here. 

Figure 3 shows that for \y\ < 1, shadowing is only a few percent effect due to the 
restricted x regions. The composition changes, with charged kaons slightly enhanced and 
protons slightly suppressed compared to charged pions. This difference is due to isospin. 
The fragmentation functions assume u and s quarks and antiquarks fragment identically to 
charged kaons [9] so that at large pt, charged kaon production is favored in AA collisions 
relative to pp. For neutral kaons, the situation is reversed. On the other hand, a n is 
twice as likely to produce a proton than a d [9] so that proton production is favored in pp 
interactions. The isospin effect is larger than the nuclear modihcations on kaons and protons 
for Pt > 7.5 GeV/c. The dominance of pion production by gluons and the assumption that 
u and d are equally likely to produce charged pions [9] leads to a negligible isospin effect 
for all Pt- Thus Raa for pions is inhuenced only by the nuclear modihcations. A similar 
conclusion was reached elsewhere [5]. 

Note that in both Figs. 2 and 3, the total Raa closely follows the pion result. This is due 
to the relative fragmentation yields: the pion yield is signihcantly higher than the proton 
yield at all pt, in contradiction with the data [1]. Reference [13] also showed that the p/n 
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ratio was underpredicted by the fragmentation functions [9]. 

Other models may predict larger shadowing effects on R\. For example, the shadowing 
model of Ref. [15], based on diffractive data from HERA, predicts large gluon shadowing for 
X < 0.01, ~ 0.3 at a; = 0.001. They also predict signihcant antishadowing for x > 0.03, 

comparable to EKS98. The gluon shadowing grows rapidly when x < 0.01. This model 
should lead to signihcant reduction in Raa for large rapidities at RHIC. However, in the 
region \y\ < 1, a; ~ 0.02 — 0.1, nuclear effects are rather small and better constrained by 
data. Thus Raa should not be signihcantly model dependent. 

Models with strong gluon saturation and/or classical gluon helds, like the colored glass 
condensate [16], predict very different nuclear gluon densities. However, px > Qs, the 
saturation scale predicted at RHIC [16]. In addition, at midrapidity, the relevant x values are 
not so small and are in a region where hxed-target data constrains the quark and antiquark 
modihcations. Theoretical analysis indicates that the gluon distributions are also rather well 
constrained in the RHIC kinematic region [17]. These studies strongly limit the possible 
effect on Raa near midrapidity. 

At the LHC, where ^/Snn = 5.5 TeV, the situation for lead on lead collisions will be 
very different. A hadron with px = 5 GeV/c (or a parton with px = 10 GeV/c) corresponds 
to a; ~ 0.002, and shadowing is signihcant, reducing Raa below 1. Away from y = 0, even 
smaller x values are probed. Gluon dominance of minijet production extends to larger px, 
considerably reducing the ehect of isospin. Thus nuclear modihcations will be signihcant at 
the LHC. 

In conclusion, nuclear shadowing cannot explain a signihcant fraction of the observed 
suppression of high-pr particles at RHIC. With the EKS98 parameterization and the nuclear 
isospin, at midrapidities 1.0 < Raa < 1-1 for charged mesons and Raa ~ 1 for protons. 
Without the restriction \y\ < 1, shadowing and isospin have a bigger ehect on Raa, with 
Raa ~ 0.7 at px = 2 GeV/c, rising to 0.8 at px = 10 GeV/c. Models with stronger 
nuclear ehects may further reduce Raa when all rapidities are considered. However, at mid¬ 
rapidity, no large ehect on Raa is expected from any model that reproduces the existing 
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lepton scattering data. 
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FIG. 1. Fractional contribution of gluons (solid curves), quarks (dot-dashed curves) and an¬ 
tiquarks (dashed curves) for (a) pion, (b) kaon, and (c) proton production in pp collisions at 
^fS = 200 GeV as a function of pT- 
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FIG. 2. Raa for pions (dashed curve), kaons (dot-dashed curve), protons (dotted curve) and 
the average over all hadrons (solid line) for gold-gold collisions at \/Smn = 200 GeV as a function 
of Pt- Spatially averaged shadowing is used. 
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FIG. 3. Raa for pions (dashed curve), kaons (dot-dashed curve), protons (dotted curve) and 
the average over all hadrons (solid line) for gold-gold collisions with |y| < 1 at \/Snn = 200 GeV 
as a function of pT- Spatially averaged shadowing is used. 
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